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Phase transitions in nuclei, small atomic clusters and self-gravitating systems de- 
mand the extension of thermo-statistics to "Small" systems. The main obstacle 
is the thermodynamic limit. It is shown how the original definition of the en- 
tropy by Boltzmann as the volume of the energy-manifold of the N-body phase 
space allows a geometrical definition of the entropy as function of the conserved 
quantities. Without invoking the thermodynamic limit the whole "zoo" of phase 
transitions and critical points/lines can be unambiguously defined. The relation to 
the Yang— Lee singularities of the grand-canonical partition sum is pointed out. It 
is shown that just phase transitions in non-extensive systems give the complete set 
of characteristic parameters of the transition including the surface tension. Nu- 
clear heavy-ion collisions are an experimental playground to explore this extension 
of thermo-statistics 

Thermodynamicists will certainly answer our question by: "Nothing" and yet: 

1 There is a lot to add to classical equilibrium statistics from our 
experience with "Small" systems: 

Following Lieb extensivity and the existence of the thermodynamic limit 
N — > oo\ N / v=const are essential conditions for conventional (canonical) ther- 
modynamics to apply. Certainly, this implies also the homogeneity of the sys- 
tem. Phase transitions are somehow foreign to this: The essence of first order 
transitions is that the systems become inhomogeneous and split into different 
phases separated by interfaces. In general are phase transitions consequently 
represented by singularities in the grand-canonical partition sum (Yang-Lee 
singularities). In the following we show that the micro-canonical ensemble 
gives much more detailed insight. 

There is a whole group of physical many-body systems called "Small" in 
the following which cannot be addressed by conventional thermo-statistics: 

• nuclei, 
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'Dividing extensive systems into larger pieces, the total energy and entropy are equal to the 
sum of those of the pieces. 
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• atomic cluster 

• polymers 

• soft matter (biological) systems 

• astrophysical systems 

• first order transitions are distinguished from continuous transitions by 
the appearance of phase-separations and interfaces with surface tension. 
If the range of the force or the thickness of the surface layers is such that 
the number of surface particles is not negligible compared to the total 
number of particles, these systems are non-extensive. 

It is common to all these examples, that the systems are non-extensive. 
For such systems the thermodynamic limit does not exist or makes no sense. 
Either the range of the forces (Coulomb, gravitation) are of the order of the 
linear dimensions of these systems, and/or they are strongly inhomogeneous 
e.g. at phase-separation. Before inventing any new type of aon-extensive ther- 
modynamics or entropy like the one introduced by Tsallisa we should realize 
that Boltzmann's definition: 



(eo a suitable small energy constant) does not invoke the thermodynamic limit, 
nor additivity, nor extensivity, nor concavity of the entropy S(E, N) (down- 
wards bending). This was largely forgotten since hundred years. We have to go 
back to pre Gibbsian times. It is a purely geometrical definition of the entropy 
and applies as well to "Small" systems. Moreover, the entropy S(E, N) as de- 
fined above is everywhere single-valued and multiple differentiable. There are 
no singularities in it. This is the most simple access to equilibrium statistics. 
We will explore its consequences in this contribution. Moreover, we will see 
that this way we get simultaneously the complete information about the three 
crucial parameters characterizing a phase transition of first order: transition 
temperature T tri latent heat per atom qi at and surface tension cr sur f. Boltz- 
mann's famous epitaph above contains everything what can be said about 
equilibrium thermodynamics in its most condensed form. W is the volume of 
the sub-manifold at sharp energy in the 6iV-dim. phase space. 



S = k * InW 



with 



W(E, N, V) = e trS(E - H N ) 



(1) 
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2 Relation of the topology of S(E, N) to the Yang-Lee singularities 

In conventional thermo-statistics phase transitions are indicated by singulari- 
ties of the grand-jCanpnical partition function Z(T, n,V), V is the volume. See 



more details in&Qtffl 
Z(T,n,V) 
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in the thermodynamic limit V — > oo\N/y =const . 

The double Laplace integral (Q) can be evaluated asymptotically for large 
V by expanding the exponent as indicated in the last line to second order in 
Ae, An around the "stationary point" e s ,n s where the linear term vanish: 

1 _ 

T ~ ~dE 
H dS 



T dN 
P__ dS_ 
T ~ 8V 



(3) 



the only term remaining to be integrated is the quadratic one. If the two 
eigen- curvatures Ai < 0, A2 < this is then a Gaussian integral and yields: 

t/2 rr 00 

Z{T,H,V) = }— e -y^-^s-Ts(e t , nts )]/T // ^ ^ ^1X^+^/2 (4) 
£0 JJo 



Z{T,p,V) = e- F(T <^ 

F(T,n,V) Tln(Vdct(e s ,n s )) InV 
> e s - nn s - Ts s + + o( — ). 

Here dct(e s , n s ) is the determinant of the curvatures of s(e, n). 

d 2 s d 2 s 



det(e, n) 



s d s 



di'Dn dn' 2 
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Ai > A 2 



(5) 
(6) 

(7) 



In the cases studied here A2 < but Ai can be positive or negative. If 
det(e s ,n s ) is positive (Ai < 0) the last two terms in eq.(||) go to 0, and we 
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Figure 1: MMMC simulation of the entropy s(e) (e in eV per atom) of a system of 1000 
sodium atoms with realistic interaction at an external pressure of 1 atm. At the energy per 
atom ei the system is in the pure liquid phase and at es in the pure gas phase, of course 
with fluctuations. The latent heat per atom is qi at = e% — e.\. Attention: the curve s(e) is 
artifically sheared by subtracting a linear function 25 + e * 11.5 in order to make the convex 
intruder visible. s(e) is always a steeply monotonic rising function. We clearly see the 
global concave (downwards bending) nature of s(e) and its convex intruder. Its depth is the 
entropy loss due to the additional correlations by the interfaces. From this one can calculate 
the surface tension cr aur f /Ttr = &s snr f * N/N sur f. The double tangent is the concave hull 
of s(e). Its derivative gives the Maxwell line in the caloric curve T(e) at Tt r 

obtain the familiar result f(T,fi,V) = e s — fin s — Ts s . I.e. the curvature Ai 
of the entropy surface s(e,n,V) decides whether the grand- canonical ensemble 
agrees with the fundamental micro ensemble in the thermodynamic limit. If 
this is the case, Z(T,fi, V) is analytical and due to Yang and Lee we have a 
single, stable phase. Or otherwise, the Yang-Lee singularities reflect anoma- 
lous points/regions of Ai > (det(e, n) < 0). This is crucial. As det(e s ,n s ) 
can be studied for finite or even small systems as well, this is the only proper 
extension of phase transitions to "Small" systems. 

3 The physical origin of the wrong (positive) curvature Ai of s(e s , n s ) 

Before we proceed with the general micro-canonical classification of the various 
types of phase transitions we will now discuss the physical origin of convex 
(upwards bending) intruders in the entropy surface. 

As nuclear matter is not accessible to us we should investigate atomic 
clusters and compare their thermodynamic behavior with that of the known 
bulk. 

In the figure ([!]) we compare the "liquid-gas" phase transition in sodium 
clusters of a few hundred atoms with that of the bulk at 1 atm.. 
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N 


200 


1000 


3000 


bulk 




T tr [K] 


940 


990 


1095 


1156 




qiat [eV] 


0.82 


0.91 


0.94 


0.923 


Na 


Sboil 


10.1 


10.7 


9.9 


9.267 






0.55 


0.56 


0.45 








39.94 


98.53 


186.6 


oo 






2.75 


5.68 


7.07 


7.41 



Table 1: Parameters of the liquid-gas transition of small sodium clusters (MMMC- 
calculation) in comparison with the bulk. 



In the table (g) we show in comparison with the known bulk values the 
four important parameters, transition temperature T tr , latent heat per atom 
qiat , the entropy gain for the evaporation of one atom Sboii as proposed by the 
empirical Trouton's rule (~ 10) and As sur f, the surface entropy per atom as 
defined above. nH^, is the average number of surface atoms of the clusters. 
a/T tr is the surface tension over the transition temperature. 



4 Phase transitions in the micro-ensemble: The topology of the 
determinant of curvatures of s(e,n): 

Now we can give a systematic and generic classification of phase transitions 
which applies also to "Small" systems and their relation to the Yang-Lee sin- 
gularities: 

• A single stable phase by d(e, n) > 
(Ai < 0). Here s(e,n) is con- 
cave (downwards bended) in both § 
directions. Then there is a one to c 
one mapping of the grand-canonical 
^->the micro-ensemble. In the two 
examples on the right the order pa- 
rameter, the direction Vi of the 
eigenvector of largest curvature Ai 
was simply assumed to be the en- 
ergy. 




enerffv 
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• A transition of first order with 
phase separation and surface tension 
is indicated by d(e, n) < (Ai > 0). 
s(e, n) has a convex intruder (up- § 
wards bended) in the direction vi c 
of the largest curvature. The whole 
convex area of {e,n} is mapped into 
a single point in the canonical en- 
semble.. I.e. if the curvature of 
S(E, N) is Ai > both ensembles 
are not equivalent. 

• A continuous ( "second order" ) transition with vanishing surface ten- 
sion, where two neighboring phases become indistinguishable. This is at 
points where the two stationary points move into one another. I.e. where 
d(e, n) — and va=o • Vd = 0. These are the catastrophes of the Laplace 
transform E — > T. Here v^ = o is the eigenvector of d belonging to the 
largest curvature eigenvalue A = (— > order parameter). 

• Finally, a multi-critical point where more than two phases become 
indistinguishable is at the branching of several lines in the {e, n}-phase- 
diagram with d = 0, Vei = 0. 

An example showing .alLpthese possible types of phase transitions in a small 
system is discussed in&oEl 
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5 On the statistical formation of multi-star systems under rotation 

Having discussed the micro-canonical description of phase transitions in "Small'' 
systems we sketch here the relevance of our new formulation of thermo-statistics 
for astrophysical systems. A finite self-gravitating system is controlled by its 
accessible phase space and its topology: 

• Due to the long-range gravity the system is non-extensive. 

• Because there is no heat-, no angular momentum bath a microcanonical 
treatment under the observation of the conservation laws is neccessary. 

• Attractive gravity leads to a collaps (phase transition of first order). 

• Under rotation the collapsed system may multi-fragment and multi-star 
systems will be formed. 
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• Also the fragmentation of Shoemaker-Levy 9 may be viewed as statistical 
multifragmcntation under the linear stress of Jupiter's gravitation field. 

With the analogy to the long-range Coulomb force one may study many aspects 
by the study of collision of heavy nuclei. E.g.: 

• At higher rotation the system prefers larger moment of inertia. This leads 
to more symmetric heaviest fragments. For the astro-problem higher 
angular momentum leads to double or multiple stars instead of a mono 
star. 

• For heavy-ion collisions, rotation leads to higher radial energy of frag- 
ments, larger variance in the statistical fragmentation which must be 
distinguished from from radial flow ! 



6 Conclusion 

Instead of using the boiling of water we used the boiling of sodium to demon- 
strate the physics of first order phase transitions and especially the surface 
tension. "Small" systems show clearly how to determine all important char- 
acteristica of first order transitions: transition temperature, latent heat, and 
surface tension. This is possible just because these systems are not extensive 
and their entropy s(e,n) has convex intruders from which the surface entropy 
can be determined. In the thermodynamic limit this fact becomes more ob- 
scured than illuminated. 

We also mentioned collapsing of self-gravitating and rotating hydrogen 
clouds towards multi-star systems. As far as the structure of the phase space 
is concerned there is a striking analogy to the phase space of fragmenting hot 
nuclei under rotation. The experimental study of the latter is also interesting 
in view of this analogy. One of the main differences is of course that the latter 
are experimentally accessible. Another important question is whether there is 
a kind of statistical equilibrium established in the astrophysical systems and 
if what the equilibration mechanism might be. Here our learning process has 
just started. 
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